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Abstract 

Background: Fragile X syndrome and tuberous sclerosis are genetic syndromes that both have a high rate of 
comorbidity with autism spectrum disorder (ASD). Several lines of evidence suggest that these two monogenic 
disorders may converge at a molecular level through the dysfunction of activity-dependent synaptic plasticity. 

Methods: To explore the characteristics of transcriptomic changes in these monogenic disorders, we profiled 
genome-wide gene expression levels in cerebellum and blood from murine models of fragile X syndrome and 
tuberous sclerosis. 

Results: Differentially expressed genes and enriched pathways were distinct for the two murine models examined, 
with the exception of immune response-related pathways. In the cerebellum of the Fmrl knockout {Fmrl-KO) 
model, the neuroactive ligand receptor interaction pathway and gene sets associated with synaptic plasticity such 
as long-term potentiation, gap junction, and axon guidance were the most significantly perturbed pathways. The 
phosphatidylinositol signaling pathway was significantly dysregulated in both cerebellum and blood of Fmrl-KO 
mice. In Fsc2 heterozygous (-1-/-) mice, immune system-related pathways, genes encoding ribosomal proteins, and 
glycolipid metabolism pathways were significantly changed in both tissues. 

Conclusions: Our data suggest that distinct molecular pathways may be involved in ASD with known but different 
genetic causes and that blood gene expression profiles of Fmrl-KO and Tsc2+/- mice mirror some, but not all, of 
the perturbed molecular pathways in the brain. 
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Background 

Autism spectrum disorder (ASD) manifests significant 
heterogeneity in part because of the interaction of under- 
lying genetic [1-3], neurobiological, and environmental 
factors [4,5] during early brain development. This hetero- 
geneity presents one of the main obstacles to the develop- 
ment of effective treatments for ASD. The complex 
genetics of ASD suggest that it is a large set of related 
disorders with diverse mechanisms; however, many of the 
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etiologies implicated in ASD may converge on a few com- 
mon pathways. Further research on single gene disorders 
associated with ASD such as tuberous sclerosis complex 
(TSC) and fragile X syndrome (FXS) may lead to an un- 
derstanding of common dysfunction at the cellular or cir- 
cuit level for a majority of ASD. In a recent survey of over 
14,000 individuals under age 35 with ASD in a Boston 
area hospital, Kohane and colleagues reported that the 
prevalence of genetic disorders of FXS and TSC in indi- 
viduals with ASD were 0.5% and 0.8% [6]. Conversely, 
30% and 50-61% of patients with FXS and TSC present 
ASD core symptoms, respectively [7,8]. If such shared 
pathophysiology exists, then treatments developed for a 
target in one disorder might be applicable to others. Mouse 
models for ASD serve an increasingly important role in 
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providing a pre-clinical test of promising pharmacological 
therapeutics [9,10]. Inactivating mutation in Tsc2 {Tsc2+I- 
mice) showed defects in axon guidance [11] and cognitive 
deficits such as impaired water maze performance [12], 
and mice with fwri -knockout (KO) presented impair- 
ments in long-term depression, hyperactivity, anxiety-like, 
and unusual social behaviors [13]. Therefore, deter- 
mining the degree to which there are shared molecular 
mechanisms in these models will inform clinical trials, 
particularly those that address populations with genetically 
heterogeneous causes of ASD. 

Although several cellular mechanisms may be impli- 
cated (reviewed in Fatemi et al. [14]), accumulating data 
support a role for the PI3K-mTOR signaling cascade in 
several genetic causes of ASD. Evidence for the PI3K- 
mTOR pathway first emerged from TSC [15,16] and 
mutations in the PTEN gene associated with ASD and 
macrocephaly [17-19]. Later, investigation of copy num- 
ber variants (CNV) in autistic individuals identified that 
PI3K-mTOR pathway-related genes were located in 
CNV hotspots [20]. These findings have led to the hy- 
pothesis that overactivation of the mTOR pathway could 
lead to abnormal synaptic function owing to an excess 
of protein synthesis at the synapse [21]. Genetic evi- 
dence that directly implicates a translation initiating 
factor, EIF4E, which is a downstream target of mTOR, 
in ASD has provided further support for this hypothesis 
[22]. Interestingly, exposure to teratogens such as val- 
proate in utero can lead to ASD in children [23], and 
valproate can also modulate this signaling pathway 
[24], suggesting that environmental factors associated 
with ASD can also play a role in PI3K-mTOR pathway 
regulation [25]. More recently, studies have found that 
PI3K-mTOR signaling is upregulated in mouse models 
of FXS, one of the most common genetic causes of 
ASD [26-28]. 

Together, the aforementioned findings suggest that an 
upregulated PI3K-mTOR signaling cascade might be a 
common mechanism in ASD and therefore would po- 
tentially be a promising drug target. Indeed, clinical 
trials using inhibitors of mTOR are already in progress 
in patients with TSC. We hypothesized that if the PI3K- 
mTOR signaling pathway is dysregulated in various 
causes of ASD, then these disorders should present with 
a similar gene expression profile signature. We chose to 
analyze TSC and FXS, two Mendelian disorders highly 
associated with ASD. Better understanding of similarities 
and differences of the cellular and molecular defects 
leading to abnormal neurological function in these two 
disorders is essential to the development of new therap- 
ies for ASD. Here, we used mouse models available for 
both genetic disorders to investigate the similarities and 
differences between gene expression profiles in the brain 
and blood cells. 



Methods 

Murine models of fragile X syndrome and tuberous 
sclerosis 

To identify molecular signatures of each mouse model 
of ASD, we performed gene expression profiling on cere- 
bellum and peripheral blood collected from two mouse 
models and compared to wild-type (WT) controls. We 
used cerebellum where the most consistent abnormalities 
were reported in the patients with ASD [14]. Post-mortem 
studies have shown a reduced number of Purkinje cells 
(PC), and several neuroimaging studies reported enlarged 
cerebella in ASD [29,30]. The cerebellum is also impli- 
cated in social interaction [31], and the loss of Tscl from 
cerebellar PC was associated with autistic-like behaviors 
[32]. Additionally, we profiled whole blood from the same 
individual mouse to compare with the gene expression 
changes in cerebellum. 

All male C57BL/6 congenic Fmrl-KO mice and 
Tsc2+/- mice with mixed 129/Sv}ae-C57BL/6 } back- 
ground have been previously described [33,34]. We pro- 
filed Tsc2+/- mice since homozygous Tsc2 KO was 
embryonic lethal. The mice were killed at 8-10 weeks 
of age following the institutional animal care and use 
committee (lACUC) euthanasia criteria (the Boston 
Children's Hospital lACUC animal protocol no. 12-07- 
2227R). For the Fmrl-KO model, 5 KO and 5 WT mice 
were profiled, and for the Tsc2+/- model 3 transgenic and 
3 WT mice were profiled. Paired blood and cerebellum 
samples were prepared for gene expression profiling. 

Genome-wide gene expression profiling using 
microarrays 

A total of 250 ng RNA was processed using established 
Affymetrix protocols for the generation of biotin-labeled 
cRNA, and the hybridization, staining, and scanning of 
arrays were performed. Briefly, total RNA was converted 
to double-stranded cDNA using a T7 primer and biotin- 
labeled cRNA was then generated from the cDNA by 
in vitro transcription. The cRNA was quantified (using 
A260) and fragmented. Fragmented cRNA was hybrid- 
ized to the Affymetrix Mouse Gene ST 1.0 array and 
scanned on an Affymetrix GeneChip scanner 3000 at 
2.5 |im resolution [35]. Microarray data are available at 
the Gene Expression Omnibus database (GSE40630). 

Validation of gene expression changes using quantitative 
RT-PCR 

Total RNA was extracted using TRIzol according to the 
manufacturer's instruction. The RNA amount was mea- 
sured using the Nanodrop (Thermo Scientific); 100 ng 
of total RNA was reversed transcribed using a cDNA 
reverse transcription kit with random primers (Applied 
Biosystems). SyBr Green PGR Master Mix (Applied Bio- 
systems) was used to amplify and detect signals from 
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cDNA with 1 mM gene-specific primers. Amplified 
signals were collected by the 7300HT Fast Real-Time 
System (Applied Biosystems) and normalized to Gapdh. 
Primer sequences used for this study were Gapdh forward 
5'-tgtgtccgtcgtggatctga-3' reverse 5'-cctgcttcaccaccttcttga- 
3', Fmrl forward 5'-ggtcaaggaatgggtcgagg-3', reverse 5'- 
agtcgtctctgtggtcagat-3', Tsc2 forward 5'-cagtgtcgac 
cagctgtctt-3', reverse 5'-tcacgctgtctggtcttgtc-3', EpsSll 
forward 5'-cagctacaacacgagaagcg-3', reverse 5'-ccgaaccttc 
caccatttgc-3 ' and GrinSa forward 5'-ctgaaacctgggtgtgaggt- 
3', and reverse 5'-aatgctgttcccacacaaca-3'. 

Microarray analysis 

All microarrays were normalized together at the probe- 
level using a quantile method, and the Affymetrix Probe 
Logarithmic Intensity ERror (PLIER) model was used to 
calculate the absolute gene expression levels as previ- 
ously described [35]. We fitted a linear model of the tissue 
(i.e., blood vs. cerebellum) and treatment (i.e., transgenic 
vs. WT) as predicting variables to each probe set. Two 
murine models were analyzed separately as different back- 
ground strains were used. Differentially expressed genes in 
each murine model were compared to the differentially 
expressed genes between the wild types of two models. 
The false discovery rate (FDR) was calculated using Storey 
and Tibshirani's method [36]. We did not use non- 
parametric tests such as the Wilcoxon rank sum test be- 
cause of the granularity of the test statistics with a small 
number of samples per group. 

We identified enriched pathways using the Gene Set 
Enrichment Analysis (GSEA) [37] . The Kyoto Encyclopedia 
of Genes and Genomes (KEGG) pathways of size 15-500 
were used for pathway analysis. Due to the relatively small 
number of samples in each group, we randomly sampled 
gene sets of equal size for each KEGG pathway to calcu- 
late the background distribution of enrichment scores. 
This procedure was done with 2,000 random drawings, 
thus the minimum permutation p-value was 0.0005. We 
used permutation /i-value 0.05 as the significance 
threshold for GSEA, and corresponding FDRs were 
described. To identify the core set of genes that ac- 
counts for the gene set's enrichment signal, we used 
leading edge analysis as described in Subramanian 
et al. [37] where a leading edge subset was defined as the 
subset of genes that gave the maximum enrichment score. 
These genes were the topmost correlated genes with 
phenotype in a gene set. 

To compare differentially expressed genes with the list 
of known ASD candidate genes as curated in the SFARI- 
genes 2.0 database (http://gene.sfari.org/) [38], we mapped 
mouse genes to human homologs using the Mouse 
Genome Informatics (MGI) Web database (http://www. 
informatics.jax.org) and performed hypergeometric tests 
to check the significance of overlap. 



Results 

Distinct gene expression changes define Fmrl and Tsc2 
transgenic models 

A total of 107 and 115 probe sets were significantly 
changed in Fmrl-KO and Tsc2+/- mice compared to 
corresponding WT littermates, respectively (uncorrected 
p-vahxe < 0.01). Not surprisingly, Fmrl was the most sig- 
nificantly downregulated gene in Fmrl-KO (p-value 
8.85 X 10'^, corresponding FDR 0.29). We used nominal 
/^-values estimated from a linear model less than 0.01 to 
rank significant probe sets because multiple testing 
correction procedures did not make any gene significant. 
The expression levels of 16 out of 107 significant probe 
sets (15.0%) in Fmrl-KO mice and 58 out of 115 signifi- 
cant probe sets (50.4%) in Tsc2+/- mice did not show a 
significant difference between blood and brain. We used 
an agglomerative hierarchical clustering with the signifi- 
cant probe sets for each model to explore the similarity 
of gene expression profiles in two genotypes and across 
tissue types. Samples were clearly separated by tissue 
type and then by genotype (Figure lA and IB). Interest- 
ingly, 71% of significant probe sets in Fmrl-KO mice 
were highly expressed in cerebellum compared to blood 
(Figure lA). For the Tsc2+/- model, the average expres- 
sion levels of 63 out of 115 probe sets (57.5%) were 
higher in cerebellum (Figure IB). Differentially expressed 
genes with statistical scores are listed in Additional 
file 1: Table 51 (Fmrl-KO vs. WT) and Additional file 1: 
Table S2 (Tsc2+/- vs. WT). 

We compared the gene expression profiles of WT 
mice since Tsc2+/- mice had 129/SvJae-C57BL/6 J back- 
grounds compared to C57BL/6 congenic backgrounds of 
Fmrl-KO mice. A total of 1,486 probes sets were differ- 
entially expressed between two WT strains. Seven probe 
sets, representing five genes, Cog7, Cc2dla, Smurfl, 
Sec31a, and AU040320, overlapped with the differen- 
tially expressed genes in Fmrl-KO vs. WT comparison. 
For the differentially expressed genes in Tsc2+/- vs. WT 
comparison, 13 probe sets (representing 9 genes: Fhadl, 
PgamS, Cts8, Piwil4, TmemlOl, Heatr7a, Mucll, 
2210404JllRik, and FamlSlb) were significantly differ- 
ent between two WT mice strains. 

Epidermal growth factor (EGF) receptor pathway 
substrate 8-like 1 (EpsSll) (Affymetrix probe set ID: 
10549655) was the only gene that was significant in 
both murine models. EpsSll was upregulated in Tsc2+/- 
mice, but downregulated in Fmrl-KO mice. Post hoc 
two-group comparison of the transgenic model to WT for 
each tissue showed that EpsSll was upregulated in blood 
(Welch's ^-test /7-value 0.07) and brain (Welch's f-test 
/?-value 0.06) of Tsc2+/- mice, while it was downregu- 
lated in blood (Welch's i-test j?-value 0.015) of Fmrl-KO 
mice. Downregulation of EpsSll in brain of Fmrl-KO was 
not significant (Welch's i-test /7-value 0.38). The homolog 



Kong et al. Molecular Autism 2014, 5:16 
http://www.molecularautism.eom/content/5/1/16 



Page 4 of 1 1 



(A) Fmrt-KO vs. wildtype (B) rsc2 +/- vs. wildtype 




Blood 
Brain 
Wildtype 
Transgenic 



Figure 1 Agglomerative hierarchical clustering of blood and brain samples. (A) Hierarchical clustering of Fmrl-KO and wild-type samples. 
One hundred seven differentially expressed transcripts (nominal p-value < 0.01) are used for hierarchical clustering of transcripts {rows in the 
heatmap) and samples [columns in the heatmap). The gene expression levels are normalized across samples, and red (upregulated) and blue 
(downregulated) are color coded according to the bottom color bar Color keys on the top of the heatmap denote tissue type and transgenic 
model. By far the most pronounced clustering is by tissue type. Within the blood samples, the two mice strains are clearly separated. (B) Hierarchical 
clustering of Tsc2 +/- and wild-type samples One hundred fifteen significant transcripts are used for cluster analysis Two mice strains formed separate 
clusters in each tissue type. The tissue specificity is not significant for the differentially expressed transcripts in Tsc2 +/- mice, whereas a majority of 
differentially expressed genes are highly expressed in cerebellum of Fmr!-KO mice (lower right cluster in the heatmap of Figure 1 A). 



of EpsSll in human, EPS8L1 is a member of EPS8-related 
proteins that play an important role in actin remodel- 
ing in response to EGF [39], and EPS8 is one of the 
reported targets of FMRP [40]. Interestingly, Stamata- 
kou and colleagues reported that Eps8 loss of func- 
tion impaired the structural and functional plasticity 
of synapses induced by long-term potentiation in pri- 
mary rat hippocampal neurons [41]. Phenotypically, 
Eps8-KO mice have impaired learning and memory, 
and excessive synaptic growth and abnormal spine 
morphology were observed in the hippocampus [42]. 
In human samples, the average expression of EPS8 in 
fusiform gyri was significantly lower among the pa- 
tients with ASD compared to controls [42]. Although 
Eps8 itself was not significantly changed in our study, 
EPS8 family genes are interesting candidates for fur- 
ther investigation. Thus, we performed qRT-PCR of 
Eps8 in brain samples of two murine models. 

Of the differentially expressed genes, several were 
also found in the expert curated database of autism 
candidate genes. Of these, Chd7, Fmrl, and Tmlhe 



were differentially expressed in Fmrl-KO mice, and 
Oxtr and Taflc were significantly changed in Tsc2+/- 
mice. The differentially expressed genes were not signifi- 
candy enriched for known ASD candidate genes in human 
(hypergeometric test p-v&lxxe 0.48 for Fmrl-KO and 0.70 
for Tsc2+/-). 

Validation of differentially expressed genes using 
quantitative RT-PCR 

Quantitative RT-PCR for individual genes was used to 
further confirm the results of the expression profiling in 
the brain samples used for the initial analysis. As expected, 
Fmrl gene expression was significantly decreased in 
Fmrl-KO mice as compared to controls (average 5.24-fold 
downregulated, Welch's i-test p-vzhxe 0.025), while Fmrl 
was unchanged in Tsc2+/- mice (Welch's i-test p-value 
0.70). Tsc2 gene expression showed a trend of decreased 
expression in both Tsc2+I- and Fmrl-KO mice as com- 
pared to WT mice that did not reach significance (average 
fold change 1.12 and 1.31 downregulation, Welch's i-test 
/?-value 0.55 and 0.27, respectively). We previously found 
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a significant reduction of Tsc2 protein in cortical neurons 
of Tsc2+/- mice [11]; however, Tsc2 mRNA expression 
was not significantly downregulated. Further, our re- 
sults confirm that EpsSll exhibited a three-fold increase in 
Tsc2+I- animals as compared to WT mice (average 3.15- 
fold upregulation, Welch's i-test /7-value 0.06). However, 
EpsSll expression shows a decreased trend — an average 
1.26 fold downregulation — in Fmrl-KO compared to WT 
similar to that observed in the expression profQing but not 
significantly changed (Welch's i-test p-vahie 0.74). 

We observed significant downregulation of subtype 3a 
of the N-methyl-D-aspartate receptor gene (GrinSa) in 
the blood of Fmrl-KO mice (uncorrected p-value 0.0098). 
GrinSa was not significant in brain with microarray data; 
however, a quantitative RT-PCR analysis of the same sam- 
ples showed a significant downregulation of this gene in 
Fmrl-KO brain (5.69-fold downregulated, Welch's f-test 
/»-value 0.046). 

Enriched biological pathways in blood and brain of the 
two mouse models 

We explored whether similar sets of biological pathways 
were perturbed in both models using GSEA [37] as there 



was only one overlapping gene — EpsSll — between the 
two lists of differentially expressed genes in Fmrl-KO 
and Tsc2+/- mice. The genes that contributed to making 
a pathway significant were identified using leading edge 
analysis. First, all genes were ranked by a per-gene 
signal-to-noise ratio that was defined as mean difference 
divided by the sum of standard deviation of each group. 
Then a running sum was calculated for each gene set. 
Beginning with the top-ranked gene, the running sum 
increased when a gene in a gene set was found and 
decreased otherwise. The enrichment score (ES) was de- 
fined to be the largest value of the running sum, and the 
genes that maximized ES were defined as the leading 
edge subset. 

Two pathways related to cytokine and complement- 
mediated signaling, cytokine-cytokine receptor inter- 
action (HSA04060, FDR < 0.0005) and complement and 
coagulation cascades (HSA04610, FDR < 0.0005), were 
the most significantly enriched pathways in Fmrl-KO 
brain (Table 1). Neuroactive ligand receptor inter- 
action (HSA04080, FDR < 0.0005), long-term potentiation 
(HSA04720, FDR 0.0023), gap junction (HSA04540, FDR 
0.0429), and axon guidance (HSA04360, FDR 0.048) were 



Table 1 Enriched pathways in blood and brain of Fmrl knockout mice 



KEGG categories 


Name 


SIZE 


NES 


NOM p-val 


FDR q-val 


Brain 












Nervous system 


HSA04080 Neuroactive ligand receptor interaction 


223 


2.34 


< 0.0005 


0.0000 




HSA04720 Long-term potentiation 


62 


-2.15 


< 0.0005 


0.0023 




HSA04540 Gap junction 


77 


-1.75 


< 0.0005 


0.0429 




HSA04360 Axon guidance 


125 


-1.81 


< 0.0005 


0.0480 


mmune system 


HSA04060 Cytokine-cytol<ine receptor interaction 


212 


2.72 


< 0.0005 


< 0.0005 




HSA04610 Complement and coagulation cascades 


47 


2.32 


< 0.0005 


< 0.0005 




HSA04650 Natural killer cell mediated cytotoxicity 


95 


1.70 


0.0015 


0.0379 




HSA04612 Antigen processing and representation 


42 


1.67 


0.0049 


0.0397 




HSA04640 Hematopoietic cell lineage 


66 


2.09 


< 0.0005 


0.0004 


Signaling patlnways 


HSA01430 Cell communication 


118 


2.14 


< 0.0005 


0.0004 




HSA04630 Jak-STAT signaling pathway 


137 


2.07 


< 0.0005 


0.0005 




HSA04070 Phospiiatidylinositol signaling system 


65 


1.77 


0.0029 


0.0420 




HSA04910 Insulin signaling pathway 


126 


1.80 


< 0.0005 


0.0455 


Metabolism 


HSA00150 Androgen and estrogen metabolism 
HSA00590 Arachidonic acid metabolism 


33 
39 


1.84 
1.81 


0.0043 
0.0017 


0.0118 
0.0147 




HSA00361 y-hexachlorocyclohexane degradation 


17 


1.69 


0.0118 


0.0370 




HSA00020 Citrate cycle 


25 


1.93 


0.0023 


0.0375 




HSA00592 Q-linolenic acid metabolism 


15 


1.70 


0.0127 


0.0392 




HSA00251 Glutamate metabolism 


30 


1.87 


0.0038 


0.0430 


Folding, sorting, and degradation 


HSA04120 Ubiquitin mediated proteolysis 


36 


1.78 


0.0039 


0.0453 




HSA04130 Snare interactions in vesicular transport 


29 


1.84 


< 0.0005 


0.0450 


Blood 












Signaling patliways 


HSA04070 Phosphatidylinositol signaling system 


65 


1.78 


< 0.0005 


0.0500 
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also enriched in Fmrl -KO brain. Among the signaling 
pathways, the PI3K signaling pathway (HSA04070) was 
changed in both brain (FDR 0.042) and blood (FDR 0.05) 
of Fmrl-KO mice. Eight genes were in the leading edges 
of the two tissues (see Methods). These were Dgkb, Dgkg, 
Dgkh, Inpp4b, InppSa, ItprS, Plcb4, and Prkca. Glutamate 
metabolism (HSA00251) was also enriched in Fmrl-KO 
brain (FDR 0.043). 

Eight pathways were significantly enriched in Tsc2+/- 
brain, while 11 pathways were enriched in Tsc2+/- 
blood (Table 2). Pathways associated with the immune 
system were significantly enriched in both blood and brain. 
Ribosome (HSA03010, FDR < 0.0005), cytokine-cytokine 
receptor interaction (HSA04060, FDR < 0.0005), and oxida- 
tive phosphorylation (HSA00190, FDR < 0.0005) were the 
most significant pathways in the brain, and the ribosome 
pathway was also significantly deregulated in Tsc2+/- 
blood (FDR 0.0468). In the Tsc2+/- brain, cytokine- 
cytokine receptor interaction (FDR < 0.0005) and 
hematopoietic cell linage (HSA04640, FDR 0.0389) 
were significant, while the Toll-like receptor signaling 
pathway (HSA04620, FDR 0.0016) and B-cell receptor 
signaling pathway (HSA04662, FDR 0.0072) were sig- 
nificant in Tsc2+/- blood. None of the differentially 
expressed genes between two WT strains was a mem- 
ber of significant pathways that we identified above. 



Common signature of Fmrl and Tsc2 transgenic models 

Cytokine-cytokine receptor interaction pathway, hematopoietic 
cell linage, and cell communication were enriched in brain 
gene expression profiles of both Fmrl-KO and Tsc2+/- 
mice (Figure 2). We performed leading edge analysis to 
find core genes that made a pathway significant, although 
each gene was not necessarily differentially expressed. For 
the cytokine- cytokine receptor pathway (N = 212), 105 
and 106 genes were the core genes in Tsc2+/- and Fmrl- 
KO brain data sets. Il7r was the only gene that showed 
marginal significances in both data sets (post-hoc Welch's 
f-test j?-values 0.011 and 0.013 in Tsc2+/- and Fmrl-KO 
brain profiles, respectively). The same gene was also the 
only common significant gene for the hematopoietic cell 
lineage pathway. Thirty-two genes were in common 
between the leading edges of cell communication for 
Fmrl-KO and Tsc2+/- brain; however, no genes were 
significantly differentially expressed. Interestingly, Tsc2 
gene expression was downregulated in Fmrl-KO blood 
(post-hoc Welch's i-test p-vahie 0.0045). 

Discussion 

We hypothesized that Fmrl-KO and Tsc2+/- mice 
would have similar gene expression profiles. In contrast, 
our findings indicate that different gene expression sig- 
natures define these two monogenic mouse models of 



Table 2 Enriched pathways in blood and brain of Tsc2+/- mice 



KEGG categories 


Name 


SIZE 


NES 


NOM p-val 


FDR q-val 


Brain 












mmune system 


HSA04060 Cytokine-cytokine receptor interaction 


212 


2.37 


< 0.0005 


< 0.0005 




HSA04640 Hematopoietic cell lineage 


66 


1.75 


< 0.0005 


0.0389 


Signaling patliways 


HSA04010 MARK signaling pathway 


241 


-1.92 


< 0.0005 


0.0431 




HSA01430 Cell communication 


118 


2.15 


< 0.0005 


0.0007 


Translation 


HSA03010 Ribosome 


66 


3.21 


< 0.0005 


< 0.0005 


Metabolism 


HSA00563 Glycosylphosphatidylinositol anchor biosynthesis 


21 


-1.85 


0.0017 


0.0444 




HSA00190 Oxidative phosphorylation 


110 


2.75 


< 0.0005 


< 0.0005 




HSA00980 Metabolism of xenobiotics by cytochrome P450 


35 


1.77 


< 0.0005 


0.0379 


Blood 












Nervous system 


HSA05010 Alzheimer's disease 


26 


1.92 


< 0.0005 


0.0034 


mmune system 


HSA04620 Toll-like receptor signaling pathway 


94 


2.00 


< 0.0005 


0.0016 




HSA04662 B cell receptor signaling pathway 


60 


1.85 


< 0.0005 


0.0072 


Cell growth and death 


HSA04110 Cell cycle 


107 


-2.31 


< 0.0005 


0.0005 




HSA04210 Apoptosis 


73 


1.78 


0.0006 


0.0181 


Translation 


HSA03010 Ribosome 


66 


1.68 


0.0025 


0.0468 


Metabolism 


HSA01032 Glycan structures degradation 


29 


2.16 


< 0.0005 


0.0000 




HSA00530 Aminosugars metabolism 


29 


1.95 


0.0007 


0.0035 




HSA00531 Glycosaminoglycan degradation 


17 


1.92 


< 0.0005 


0.0043 




HSA0051 1 Other glycan degradation 


15 


1.88 


0.0007 


0.0047 




HSA00600 Sphingolipid metabolism 


31 


1.69 


0.0064 


0.0473 



Kong et al. Molecular Autism 2014, 5:16 
http://www.molecularautism.eom/content/5/1/16 



Page 7 of 1 1 



KEGG categories 
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HSA04060 Cytokine-cytol<ine receptor interaction 
HSA04640 Hematopoietic cell lineage 
HSA04650 Natural killer cell mediated cytotoxicity 
HSA0461 2 Antigen processing and representation 
HSA04610 Complement and coagulation cascades 
HSA04620 Toll-like receptor signaling pathway 
HSA04662 B cell receptor signaling pathway 

HSA04080 Neuroactive ligand receptor interaction 
HSA04720 Long-term potentiation 
HSA04360 Axon guidance 
HSA04540 Gap junction 
HSA05010 Alzheimer's disease 



■ 
1 



Translation 



HSA03010 Ribosome 



Signaling pathways 



Cell growth and death 



HSA01430 Cell communication 
HSA04070 Phosphatidylinositol signaling system 
HSA04630 Jak-STAT signaling pathway 
HSA04910 Insulin signaling pathway 
HSA04010 MARK signaling pathway 

HSA04110 Cell cycle 
HSA04210 Apoptosis 



Folding, sorting, and degradation HSA04120 Ubiquitin mediated proteolysis 

HSA04130 Snare interactions in vesicular transport 



IVIetabolism 



HSA00020 Citrate cycle 

HSA00361 Y-hexachlorocyclohexane degradation 

HSA00592 a-linolenic acid metabolism 

HSA00590 Arachidonic acid metabolism 

HSA00150 Androgen and estrogen metabolism 

HSA00251 Glutamate metabolism 

HSA0051 1 Other glycan degradation 

HSA00530 Aminosugars metabolism 

HSA00531 Glycosaminoglycan degradation 

HSA00563 Glycosylphosphatidylinositol anchor biosynthesis 

HSA00600 Sphingolipid metabolism 

HSA01032 Glycan structures degradation 

HSA00980 Metabolism of xenobiotics by cytochrome P450 

HSA00190 Oxidative phosphorylation 



Figure 2 Enriched pathways in Fmrl-KO and Tsc2+/- models of ASD, found using the Gene Set Enrichment Analysis (GSEA). Red 

(upregulated) and blue (downregulated) squares in tlie matrix represent enriched pathways for each data set (false discovery rate < 0.05). Two 
immune system pathways (cytokine-cytokine receptor signaling pathway and hematopoietic cell lineage) and one signaling pathway 
(cell communication) were significant in the brain gene expression profiles of both mice models. 



ASD. Global expression profiles of the two models ex- 
amined were distinct, such that only one gene — EpsSll — 
was in common. This is a particularly surprising result 
since translational dysregulation as well as aberrant syn- 
aptic protein synthesis associated with both disorders 
has been proposed as one possible pathway leading to 
autistic pheno types, including cognitive impairment [21]. 



Our data indicate that FXS and TSC may have very 
distinct brain and blood cellular phenotypes despite the 
fact that both syndromes result in similar behavioral and 
cognitive symptoms. 

Nonetheless, we did find that the cytokine and com- 
plement signaling pathways are differentially regulated 
in both mouse models although the specific genes 
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affected within these pathways were different. The 
immune system has been implicated in ASD in multiple 
ways, but the exact mechanism of the interaction 
between the immune system and genetic disorders that 
result in an increased risk of autism has not been well 
studied. Modulation of the immune system may not be 
completely unexpected since the TSC-mTORCl path- 
way regulates inflammatory responses after bacterial 
stimulation in monocytes, macrophages, and primary 
dendritic cells [43], and this pathway contributes to 
cytokine upregulation in response to endotoxins [44]. 
However, this is of particular interest in light of a recent 
study that showed that immune activation during gesta- 
tion can markedly worsen the neurological phenotype of 
Tsc2+/- mouse pups [45]. Similarly, preliminary obser- 
vations indicate that plasma protein levels of a number 
of cytokines differ between individuals with and without 
FXS. Furthermore, differences in cytokine and other 
immune-signaling genes were observed between the FXS 
group with autism and the FXS group without autism 
[46]. On the other hand, Yuskaitis and colleagues inves- 
tigated the peripheral immune system of Fmrl -KO mice, 
but did not find any differences in either the T-cell popu- 
lation at basal and stimulated status or the proinflamma- 
tory cytokines TNFa and IFNy at basal and stimulated 
status [47]. How the loss of FMRP leads to changes in the 
immune system, however, remains unclear. 

Accumulating evidence over the last few years indi- 
cates that the TSC and FMRP pathways interact. How- 
ever, precisely how these pathways interact remains an 
open question. On the one hand, FMRP can be phos- 
phorylated by S6K, an enzyme downstream of TSC [48] . 
On the other hand, mTOR, the kinase inhibited by the 
TSC2 protein, has increased activity in Fmrl -KO neu- 
rons [26], and FMRP-deficient cells display increased 
activity of PI3K, an enzyme upstream of TSC proteins 
[27]. Recently Auerbach et al. reported that the synaptic 
dysfunction in the CAl region of the hippocampus of 
Tsc2+/- mice was opposite to that of Fmrl-KO mice [49]. 
In fact, manipulating the mGluR receptors with positive 
allosteric modulators was sufficient to rescue this defect in 
Tsc2+/-, while inhibiting the mGluR receptors was neces- 
sary in the case of Fmrl-KO mice. Finally, a genetic cross 
of the Tsc2+/- and Fmrl-KO mice was similar to WT in 
CAl synaptic physiology and contextual learning. These 
results argue that loss of Tsc2 and Fmrl have some oppos- 
ite cellular phenotypes, which can be rescued in a double 
knockout. However, the mechanisms by which Tsc2 and 
Fmrl result in opposite synaptic phenotypes is not yet 
clear. In post-hoc analysis, we found that Tsc2 gene ex- 
pression was downregulated in Fmrl-KO blood (Welch's 
f-test ;7-value 0.0045). Thus, in the Tsc2+/- Fmrl-KO 
mice, the Tsc2 expression level may be closer to WT and 
may contribute to the rescue of the synaptic physiology. 



This finding points to another level of interaction be- 
tween the FMRP and Tsc2 functions in the cell. Fu- 
ture studies are required to understand whether Fmrl 
loss leads to changes in Tsc2 mRNA via transcriptional or 
post-transcriptional regulation. 

Similarities between the proposed roles of TSCl/2 and 
FMRP proteins in regulation of protein synthesis have 
led to the popular hypothesis that hyperactive mTOR 
signaling is pathogenic in both FXS and TSC. Our re- 
sults indicate the gene expression dysregulation differs 
markedly in the two conditions. This is consistent with 
the previously reported changes in neuronal morphology 
in each of these models. Neurons deficient in Tscl or 
Tsc2 display lower dendritic spine density in contrast to 
Fmrl-KO neurons [50], which have increased spine 
density [51]. At the biochemical level there are also some 
important differences. In rSCi/2-null cells, mTORCl- 
dependent negative feedback mechanisms exist to dampen 
the activation of upstream components of the network 
such that Akt activation is decreased [52]. However, in 
Fmrl-KO neurons, Akt activity is enhanced [27]. While 
mTORCl may be activated because of loss of either Tscl/2 
or FMRP, the neuronal phenotype and gene expression 
profiles may be altered by changes in the activation of 
other signaling pathways. Such differences have implica- 
tions for targeted treatment options of the two distinct 
genetic conditions. 

The mouse models we have explored have been devel- 
oped as models of ASD with known divergent genetic 
etiologies. The overall divergence in gene expression 
dysregulation in these two models does not rule out 
shared downstream effects, and indeed we observed an 
overlapping dysregulation of the cytokine signaling path- 
way. Nonetheless, it does suggest that a multiplicity of 
therapeutics will have to be developed for the varied 
mechanisms contributing to the increasingly fine-grained 
distinctions between the etiologies of ASD. 

We could identify similar sets of biological pathways 
enriched in both tissues. In the Fmrl-KO mice, the PI3K 
signaling pathway was dysregulated in both blood and 
brain, while the ribosome pathway was dysregulated in 
both tissues of Tsc2 +/-. With the two mouse models, we 
also could identify biological pathways that were positively 
correlated with genetic background in both tissues. Previ- 
ous studies demonstrated that peripheral blood expression 
signatures could be used to classify the clinical conditions 
of brain disorders [35,53]. Likewise, our results suggest 
that peripheral blood signatures could be used to identify 
genotypes as well as some transcriptional changes present 
in brain. 

The current study is limited by the different back- 
ground strains of Fmrl-KO and Tsc2+/- mice and by 
the small sample size. The differentially expressed genes 
between two WT background strains overlapped with the 
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significant genes in each murine model - 5 for Fmrl-KO 
and 9 for Tsc2+/-. Although 1,486 probe sets were signifi- 
cantly differentially expressed between WT mice of two 
different backgrounds, none of these genes was a member 
of the significant pathways that characterized Fmrl-KO 
and Tsc2+/- mice in blood and brain. Due to the small 
sample size, we did not have enough statistical power for 
microarray experiments to detect and compare relatively 
low-expressed genes. Further study using the same back- 
ground strain and increasing the sample size will be essen- 
tial to confirm our findings, and using a more sensitive 
quantification method such as RNA-seq will improve sen- 
sitivity for low abundance transcripts. For the Fmrl-KO 
model, we did not include female mice heterozygous for 
Fmrl in this experiment. Heterozygous Fmrl female 
mice should exhibit genetic mosaicism due to random 
X-inactivation of one X chromosome during develop- 
ment. For this reason, most previous studies characterized 
male Fmrl KO models. Qin and colleagues performed an 
interesting comparison of male Fmrl-KO and homozy- 
gous and heterozygous KO in female mice [54]. They 
reported that only the homozygous mice had a deficit on 
the passive avoidance test, whereas both homozygous and 
heterozygous female mice exhibited hyperactivity and 
increased susceptibility to seizures. A follow-up experi- 
ment with both sexes and different dosages of Fmrl in 
female mice with a larger sample size would be ideal 
since a gender effect on global gene expression profiles 
should be considered when both sexes are included in the 
experiment. 

Conclusions 

Contrary to our initial hypothesis that Fmrl-KO and 
Tsc2+/- mice would share a transcriptional signature, 
we found that the two mouse models presented distinct 
sets of differentially expressed genes. In retrospect, this 
is not surprising as multiple lines of evidence suggest 
that FXS and TSC are actually driven by opposite mo- 
lecular phenotypes [49]. Despite these gene-level differ- 
ences, however, we observed that cytokine signaling, cell 
communication, and hematopoietic cell lineage genes 
were differentially expressed in both mouse strains. Sec- 
ond, our results show that blood expression signatures 
mirror many aspects of the brain transcriptome. Specif- 
ically, several pathways were dysregulated in both the 
brain and blood of the two mouse models studied here. 
This confirmation is important for the future use of 
blood tissue to study neurodevelopmental disorders. 

Availability of supporting data 

The data set supporting the results of this article is avail- 
able in the Gene Expression Omnibus repository with 
the accession identifier GSE40630 (http://www.ncbi.nlm. 
nih.gov/geo/query/acc.cgi?acc=GSE40630). 
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